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ABSTRACT 

A numerical study of QED radiative corrections for the SMC experiment 
at SPS is been performed. The semi-analytical program POLRAD is been used 
to get the size of the radiative corrections for the proton to get the size of the 
radiative corrections for the measurements of the proton spin dependent struc- 
ture functions gi{x) and g2{x). A brief description of the program POLRAD 
is given. 

1. Introduction 

The asymmetry of deep inelastic scattering of polarized leptons on polarized 
nucleons 

l{h)+N{pi)-^l{k2)+X{p2), (1.1) 

gives important information about the spin dependent structure functions of the 
nucleons gi^x^Q"^) and g2{x,Q'^) The experimental study of the asymmetry by 
the Spin Muon Collaboration (SMC) is based on the measurements of the scattering 
of polarized muons on polarized deuterons and protons. In the data analysis of 
this experiment the Fortran program POLRAD^ is used for the calculation of the 
radiative corrections (RC) from higher order QED processes (Fig. 1) for the cross 
section asymmetries 

^11 = ^-lu^ (1-2) 

and 

A± = (1.3) 

where a^^, a'^^ ( cr^^, cr^^) are the cross sections for inclusive deep inelastic scattering 
of longitudinally polarized muons off longitudinally (transversely) polarized nucleons. 



The radiative corrections are important for the higher statistics experiments and are 
dependent of the methods used to extract the wanted cross section from the data. 

The fixed target experiments on deep inelastic lepton-nucleon scattering mainly 
use the energy E' and the angle 9 of the scattered lepton to determine the kinematics 
of the events, i.e. the analysis of deep inelastic events is based on the measurement 
of the familiar leptonic variables 

gf = -(k, - k2)\ yi = pi(h - k2)/pih, XI = Qll{Syi), (1.4) 

with 

S ^{ki+pif ^2ME. (1.5) 

Here E is the energy of the incident lepton in the laboratory system, pi=0, and M 
is the mass of the target. 

The methods for the calculation of the RC for deep inelastic scattering of leptons 
on unpolarized nucleons in terms of the leptonic variables are well known. Both 
methods are based mainly on model-independent approaches. The model-independent 
treatment of the RC uses the phenomcnological structure functions of nucleons, which 
are defined in the Born approximation. A detailed comparison^ has shown that the 
agreement between the two methods is better than 2 %. 

The model-independent radiative corrections in terms of the variables xi, Qj to 
the cross section of deep inelastic scattering of polarized leptons by longitudinally 
polarized nucleons were calculated in.^ The formulae derived in these papers were 
generalized in^*' to take into account the case of a transversely polarized target and 
the scattering on polarized spin-1 targets. 

A new development for the model-independent approach in the calculations of 
the RC for deep inelastic scattering has been arised for the experiments at HERA 
where for the physical analysis of deep inelastic events one could use not only the 
leptonic variables (1.4) but also the kinematical variables derived from the hadron 
measurements 

Qh = -(P2 - Pl)^ Vh = Pi(p2 - Pi)/Piki, Xh = Ql/{Syh), (1.6) 
or some mixture of both. 

The theoretical investigations of the radiative corrections for the experiments 
at HERA have been summarized in the Proceedings of the Workshop "Physics at 
HERA",^^ and the different programs for the calculation of the RC have been cross- 
checked.^^ 

The model-independent treatment of leptonic QED corrections for neutral 
current ep-scattering, including the bremsstrahlung process 

l{k,) + N{p^) ^ l{k2) + X{p2) + ^{k), (1.7) 

the elastic radiative tail 



l{k,) + N{p^) ^ l{k2) + N{p2) + 7(A:). 



(1.8) 



and QED vertex corrections have been done in terms of several variables - leptonic, 
hadronic, mixed and using the Jacquet-Blondcl variables. 

The new formulae were implemented in the program package TERADQl.^^ The 
program TERAD91 can be applied to calculate the RC factor of the order 0{a) 
to the measured differential cross section, cPa^^"'^ / dxdQ"^ . After the acceptance and 
smearing corrections it reduces to 

dFa"/dxdQ'' = d^a^'^''ldxdQ^{\ + Q^)). (1.9) 



The semi-analytical program TERAD91 has been used ^^'^^ to evaluate the ra- 
diative corrections for the first measurement of the proton structure function F^ix^ Q^) 
in the range x = 10~^ — 10~^ and 5 GeV"^ < < 80 GeV"^ using the electron and 
the mixed kinematical variables }^ 

The radiative correction factor 5{x,Q'^) for the inclusive cross section of deep 
inelastic scattering is defined by the equation 



Six.Q'' 



- 1, 



:i.io) 



where d^ a^^^°^ / dxdQ"^ is the theoretical approximation for the measured cross section, 
^2^mcasy'^^^Q2^ which coutaius contributions from higher order QED and electroweak 
processes, and ^ a^°^^ j dxdQ"^ is the Born cross section of the process (1.1). 

In this note we describe the semi- analytical program POLRAD which has been 
used to calculate the radiative corrections for the measurements of the spin dependent 
proton structure functions g\{x) and gy^x^ by the SMC experiment. 

2. The semi-analytical program POLRAD 

For <^ Ml only electromagnetic processes are relevant, and the unpolarized 
part of the Born cross section of the deep inelastic scattering (1.1) is determined by 
two structure functions F^ix^Q'^^ and 



2xF^{x,Q) 



2, vl + p')i^2(a:,g2^ 



1 + i?(a;,g2) ' 



dxdQ^ g^x 



\-y + 



Here and below 
and 



X = Xi, 



2(i + i?(x,g2))_ 
y = yi, = Qf 

g2 2Mx 



Ey 



(2.1) 
(2.2) 
(2.3) 

(2.4) 



The spin dependent part of the Born cross section has contributions from the 
structure functions gi and g2 as: 



-y 



when the beam and target spins are coUinear, and 



(2.5) 



-yp\ 



y 



p2y2 



^9i{x) + g2{x) 



(2.6) 



when the spin directions of beam and target are orthogonal. 

The program POLRAD contains the formulae for leptonic QED corrections to 
the cross sections (2. 2), (2. 5) and (2.6). These corrections include photonic brems- 
strahlung from leptons, vertex corrections (see Fig.l, only 7-exchange ) and vacuum 
polarization of the order 0{q) with soft-photon exponentiation. This set of the higher 
order Fcynman graphs is standard for the model-independent treatment of the QED 
RC to deep inelastic scattering^ and gives the numerically largest contribution. 

The general formula for the radiative cross sections expressed in terms of lep- 
tonic variables can be written as a sum of two parts: 



j2 0{a) j2 Born ^ ^ J2 Born 

W ^ W <^ ^ .™(.,e^))+//d...«^(.,Q^..«) 1^ (2.7) 

The first part of (2.7), proportional to S^^, contains the contributions from vertex 
corrections, vacuum polarization and from the soft part of the real photon radiation. 
The second part accounts for the bremsstrahlung of hard photons. It depends on the 
structure functions, not only in a given F{x, Q^) point , but in the range of {xh, Ql) 
defined by kinematics of Fig. 2a. 

For radiative events 

X < Xh < I, < Ml < W'^, (2.8) 

were the square of the invariant mass of the hadronic final state can be defined 

as 

Ml -M^ + Ql (i^) , (2.9) 

\ Xfi / 

and the invariant mass 

W^^M^ + Q^(^^), (2.10) 



X 

corresponds to the hadronic jet without radiation. 

From the experience of the studies of the RC for deep inelastic scattering'^'^^ it 
is known that the large radiative corrections are mainly due to the emission of hard 



photons ( the second part of (2.7)). For hard photon bremsstrahlung in the process 
(1.7) and (1.8) the kinematic minimum of Qf^ is 

{QDmin - X^M\ (2.11) 

This formula shows that it is impossible to calculate the RC without the knowledge 
of the nucleon structure functions in the region — > 0. 

The program POLRAD, instead using of the usual M| and Q|, or Xh and 
(Fig. 2a) is using other invariant kinematical variables (Fig. 2b): 

R^W^-Ml + Q^-Ql r= ^,_g;^l_^, . (2.12) 

for the calculation of the two-dimensional integral in (2.7)- the contribution of the 
inelastic radiative tail. The boundaries of these variables are: 

< R < — , (2.13) 



1 + T 



and 



' max.rmn 



2M2' V Sy 

The point F{x, Q"^) of the Fig.2a in which we want to calculate the RC has been 
transformed into the line R = from Tmm to Tmax of the Fig. 2b, and the Born 
kinematics of the process (1.1) is defined on this line. The contribution of the elastic 
radiative tail (1.8) is determined by the integral on the line it! = {W'^ — M^) /(l + r). 

Now we briefly explain how the program POLRAD works. For any point F{x, Q"^) 
of deep inelastic scattering the POLRAD routine CONKIN evaluates the Born kine- 
matics (2.12)-(2.14) and routine BORNIN gives the Born cross sections (2.2)-(2.6). 
DELTAS calculates the factorized part of the RC 5^^ which is independent on the 
polarization of the beam or target. The two-dimensional integral in (2.7) is calculated 
in terms of the variables (2.12) by the routines QQT, QQINT, TAILS and FFU. The 
routines STRF and COMFST contain the parameterizations of the structure func- 
tions Fi^2 , 91,2 and also elastic and quasi-elastic form factors for the calculation of 
their contribution to the process (1.8). In the last two routines there is the possibility 
of putting other parameterizations of structure functions and form factors. The main 
routine POLRAD calculates the final results - the RC to the asymmetries. 



3. Numerical results and conclusion 



In this section we will give the results of the calculations of the RC with the 
help of the program POLRAD to the longitudinal A| and transverse A'l^ asymmetries 
in inclusive deep inelastic scattering of polarized muons off polarized protons. 

In order to calculate the RC one has to use the realistic parameterizations of 
the structure functions F12 and gi^2 over the full range of x and Q^. Several different 



parameterizations of the unpolarized structure functions, F^a;, Q^), including the 
NA37 parametrization have been used. For the region < 0.2 GeV^ the 
parameterizations of F2{x,Q'^) were taken at — 0.2 GeV'^ and multiphed by the 
suppression factor^" (1 — exp^—aQ"^)) , where a = 3.37 GeV^^. The Ff(a;,(5^) is 
calculated by (2.1), where the values of R{x, Q"^) were taken from the fit of the SLAC 
data .^^ 

The longitudinal spin structure function, can be obtained from the asym- 
metry A\ by the relation 

31^"^^^ >- 2x[l + R{x,Q^)]' ^'^■^> 

where Afl is the virtual photon asymmetry which can be measured experimentally. 
But here we will use the model parameterizations for from.^^"^"^ 

For the calculation of the transverse asymmetry, g2 is needed. It can be 
estimated by equation (23) of ref.^ if we neglect the contribution of the last term of 
this equation. 

Figures 3-9 shows the results for the quantities 

where AA^ is the difference between the asymmetries with and without the radiative 
corrections 

AAP = AP{Born + QED RC) - AP{Born). (3.3) 
and D\\ ± are the depolarization factors 



" y2 + 2(l-y)(l + i?)' ^ ^2 + 2(l-y)(l + i?)- 

The calculations presented here have been performed with the NA37 parametriza- 
tion^s for F2{x,Q'^). 

From Figures 5 and 6 it is seen that for the longitudinal polarization of the 
protons wc get different results for Ay4j|/Z)|| and we see a strong dependence of 
the RC on the choice of the Qi parametrization. The RC are not small at low x : 
I A/ll/Z^ii I < 0.025. At low X the influence of the hard bremsstrahlung from the 
processes (1.7) and (1.8) on the value of the RC is very big. For the results which 
were obtained from POLRAD in the case of the transverse polarization of the target 
we have used the Qi parametrization from^^ and have calculated g2 from eq. 23 of.^ 
The transverse asymmetry with such input for gi 2 in the Born approximation looks 
good for the who kinematical region. But for the QED corrected asymmetry ^4^/ D^^ 
(Fig. 8) and for absolute correction AAF^/D^^ (Fig. 9) we get results which cannot be 
understood. 



It is worth mentioning that the treatment of the radiative corrections should 
be crosscheck with the results of other programs. From this point of view the com- 
parison of the results of the semi-analytical program POLRAD with the results of 
other existing program is needed. For instance, with the Monte Carlo event generator 
HERACLES^^ where the independent calculations have been implemented.^^ 
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Fig. 1. Feynman graphs contributing in the higher order to deep inelastic IN- 
scattering. 




Fig. 2. Integration region of {xh, QD in (2-7) and the correspondent region in the 

terms of (r, R) . 




Fig. 3. Longitudinal asymmetry Ay D\\ in the Born approximation for beam energy 

190 Gey 2; gl from 




Fig. 7. Transverse asymmetry A^/D^ in the Born approximation for beam energy 

100 GeV^; g{ from .^^ 



